1 



The chromo-Weibel instability 



O 
O 

(N 

o 

Q 

> 



o 

Oh 
I 

Oh 

CD 



X 



Michael Strickland 

Frankfurt Institute for Advanced Studies 
Johann Wolfgang Goethe - Universitdt Frankfurt 
Max-von-Laue-Strafie 1 
D-60438 Frankfurt am Main 
Germany 

I discuss the physics of non-Abelian plasmas which are locally anisotropic in momentum space. Such 
momentum-space anisotropies are generated by the rapid longitudinal expansion of the matter created in the first 
1 fm/c of an ultrarelativistic heavy ion collision. In contrast to locally isotropic plasmas anisotropic plasmas 
have a spectrum of soft unstable modes which are characterized by exponential growth of transverse chromo- 
magneticZ-electric fields at short times. This instability is the QCD analogue of the Weibel instability of QED. 
Parametrically the chromo-Weibel instability provides the fastest method for generation of soft background 
fields and dominates the short-time dynamics of the system. The existence of the chromo-Weibel instability 
has been proven using diagrammatic methods, transport theory, and numerical solution of classical Yang-Mills 
fields. I review the results obtained from each of these methods and discuss the numerical techniques which are 
being used to determine the late-time behavior of plasmas subject to a chromo-Weibel instability. 
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I. INTRODUCTION 

With the ongoing ultrarelativistic heavy-ion collision ex- 
periments at the Relativistic Heavy-Ion Collider (RHIC) and 
planned the Large Hadron Collider (LHC) physicists hope 
to produce and study the properties of a thermalized quark- 
gluon plasma (QGP) which is expected to be formed when 
the temperature of nuclear matter is raised above its criti- 
cal value, T,. - 200 MeV - lO'^ K. Given the small size 
and short lifetime of the matter created in an ultrarelativistic 
heavy-ion collision this is not trivially accomplished. One of 
the chief obstacles to thermalization in ultrarelativistic heavy- 
ion collisions is the rapid longitudinal expansion of the matter 
created in the central rapidity region. If the matter expands 
too quickly then there will not be sufficient time for its con- 
stituents to interact and thermalize. During the first 1 fm/c 
after the nuclear impact the longitudinal expansion causes the 
created matter to become much colder in the longitudinal di- 
rection than in the transverse directions [1], corresponding to 
(Pl) ^ (Pt) local rest frame. After this initial period 

of longitudinal cooling, the expansion slows and one can then 
ask what are the dominant mechanisms for driving the system 
towards an isotropic thermal QGR 

The evolution of the partonic matter created during a high- 
energy nuclear collision was among the questions which the 
"bottom-up" thermalization scenario [1] attempted to answer. 
For the first time, it addressed the dynamics of soft modes 
(fields) with momenta much below Qs coupled to the hard 
modes (particles) with momenta on the order of Qs and above 
[2-4]. However, it has emerged recently that one of the 
assumptions made in this model was not correct. The de- 
bate centers around the first stage of the bottom-up scenario 
(1 <C Qil ^ a^^/^) in which it was assumed that (a) col- 
lisions between the high-momentum (or hard) modes were 
the driving force behind isotropization and that (b) the low- 
momentum (or soft) fields acted only to screen the electric 
interaction. In doing so, the bottom-up scenario implicitly as- 



sumed that the underlying soft gauge modes behaved the same 
in an anisotropic plasma as in an isotropic one. 

However, to be self-consistent one must determine the col- 
lective modes which are relevant for an anisotropic plasma 
and use those. In the case of gauge theories this turns out 
to be a qualitative rather than quantitative correction since in 
anisotropic QCD plasmas the most important collective mode 
corresponds to an instability to transverse chromo-magnetic 
field fluctuations [5-8]. This instability is the QCD analogue 
of the QED Weibel instability [5]. Recent works have shown 
that the presence of the chromo-Weibel instability is generic 
for distributions which possess a momentum-space anisotropy 
[9-1 1] and have obtained the full non-Abelian hard-loop (HL) 
action in the presence of an anisotropy [12]. Another impor- 
tant development has been the demonstration that the chromo- 
Weibel instability also exists in solutions to pure classical 
Yang-Mills fields in an expanding geometry [13-15]. 

Recently there have been significant advances in the un- 
derstanding of non-Abelian soft-field dynamics in anisotropic 
plasmas within the HL framework [16-19]. The HL frame- 
work is equivalent to the collisionless Vlasov theory of eikon- 
alized hard particles, i.e. the particle trajectories are assumed 
to be unaffected (up to small-angle scatterings with 9 ~ ^) 
by the induced background field. It is strictly applicable only 
when there is a large scale separation between the soft and 
hard momentum scales. Even with these simplifying assump- 
tions, HL dynamics for non-Abelian theories is complicated 
by the presence of non-linear gauge-field interactions. 

These non-linear interactions become important when the 
vector potential amplitude is on the order of Anon-Abeiian ~ 
Ps/g ^ fhPh, where pi, is the characteristic momentum of the 
hard particles, e.g. pf, ^ Qs for CGC initial conditions, //, is 
the angle -averaged occupancy at the hard scale, and ps is the 
characteristic soft momentum of the fields (ps ~ gfhPh)- In 
QED there is no such complication and the fields grow expo- 
nentially until AAbeiian ^ Ph/ g which point the hard parti- 
cles undergo large-angle deflections by the soft background 
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field which rapidly isotropize the system. In fact, in QED 
the Weibel instability is the fastest process driving plasma 
isotropization. In QCD, however, the effect of the non-linear 
gauge self-interactions affects the system's dynamics primar- 
ily slowing down instability-driven particle isotropization. 

To include the effects of gauge self-interactions numeri- 
cal studies of the time evolution of the gauge-covariant HL 
equations of motion are required. Recent numerical studies of 
HL gauge dynamics for SU(2) gauge theory indicate that for 
moderate anisotropics the gauge field dynamics changes from 
exponential field growth indicative of a conventional Abelian 
plasma instability to linear growth when the vector potential 
amplitude reaches the non-Abelian scale, Anon-Abeiian ~ P\i 
[17, 18]. This linear growth regime is characterized by a 
turbulent cascade of the energy pumped into the soft modes 
by the instability to higher-momentum plasmon-like modes 
[20, 21]. These results indicate that there is a fundamental 
difference between Abelian and non-AbeUan plasma instabil- 
ities in the HL limit. 

In addition to numerical studies in the HL limit there have 
been numerical results from the solution to the full non-hnear 
Vlasov equations for anisotropic plasmas [22, 23]. This ap- 
proach can be shown to reproduce the HL effective action in 
the weak-field approximation [24—26] ; however, when solved 
fully the approach goes beyond the HL approximation since 
the full classical transport theory also reproduces some higher 
n-point vertices of the dimensionally reduced effective action 
for static gluons [27]. Numerical solution of the 3d Vlasov 
equations show that chromo-instabilities persist beyond the 
HL limit [23]. Furthermore, the soft field spectrum obtained 
from full Vlasov simulations shows a cascade or "avalanche" 
of energy deposited in the soft unstable modes in higher mo- 
mentum modes similar to HL dynamics. 

II. ANISOTROPIC GLUON POLARIZATION TENSOR 

In this section I consider a quark-gluon plasma with a par- 
ton distribution function which is decomposed as [9] 

/(p) = 2Ar^[n,(p)+«^(p)] +4iV,n,(p) , (1) 

where tiq, rig, and ng are the distribution functions of quarks, 
anti-quarks, and gluons, respectively, and the numerical co- 
efficients collect all appropriate symmetry factors. Using the 
result of Ref. [9] the spacelike components of the HL gluon 
self-energy for gluons with soft momentum (k ^ g phard) can 
be written as 

(2) 

where K = (co,k), V — (l,p//?), and the parton distribution 
function /(p) is, up to an integrability requirement, com- 
pletely arbitrary. In what follows we will assume that /(p) 
can be obtained from an isotropic distribution function by the 
rescaling of only one direction in momentum space. 

In practice this means that, given any isotropic parton distri- 
bution function fiso{p), we can construct an anisotropic ver- 
sion by changing the argument of the isotropic distribution 



function, /(p) = /iso (^\/p^ + ^(p-n)^j , where n is the direc- 
tion of the anisotropy, and ^ > — 1 is an adjustable anisotropy 
parameter with ^ = corresponding to the isotropic case. 
Here we will concentrate on ^ > which corresponds to a 
contraction of the distribution along the n direction since this 
is the configuration relevant for heavy-ion collisions at early 
times, namely two hot transverse directions and one cold lon- 
gitudinal direction. 

Making a change of variables in (2) it is possible to inte- 
grate out the |p|-dependence giving [9] 

x^S'-' + ^^^V (3) 
V K-V + ieJ' ^ ' 

where cos 6„ = k • n and m|, > 0. The isotropic Debye mass, 
niD, depends on /iso but is parametrically mo ~ gphard- 

The next task is to construct a tensor basis for the space- 
Uke components of the gluon self-energy and propagator. We 
therefore need a basis for synmietric 3-tensors which depend 
on a fixed anisotropy 3-vector n' with — \. This can 
be achieved with the following four component tensor ba- 
sis: A'i = 5'j - W/l^, B'j = m/k^, C'j = ffhi/n^, and 
D'^ = k'fij + k'fi' with «' = A'^hK Using this basis we can 
decompose the self-energy into four structure functions a, p, 
Y, and 5 as 11 = a A -|- PB -|- yC -I- 6D. Integral expressions for 
a, P, Y, and S can be found in Ref. [9] and [11]. 



ni. COLLECTIVE MODES 

As shown in Ref. [9] this tensor basis allows us to express 
the propagator in terms of the following three functions 

where 2Q.I = Q? ± {O.'^ -4{{a + j+ k^)^ - k^h^b^)) ^'^ and 

^2 = a-FP-|-Y+A:2. 

Taking the static limit of these three propagators we find 
that there are three mass scales: m± and ma- In the isotropic 
limit, ^ — > 0, = ~ and = m^. However, for ^ > 
we find that < for all |9„ 1 7^ 7t/2 and < for all 
|9„ I < 7l/4. Note also that for ^ > both and have 
there largest negative values at 9„ = where they are equal. 

The fact that for ^ > both and m?. can be negative in- 
dicates that the system is unstable to both magnetic and elec- 
tric fluctuations with the fastest growing modes focused along 
the beam line (9„ = 0). In fact it can be shown that there are 
two purely imaginary solutions to each of the dispersion re- 
lations Ajj' (/T) = and Al' {K) = with the solutions in the 
upper half plane corresponding to unstable modes. We can 
determine the growth rate for these unstable modes by taking 
CO ir and then solving the resulting dispersion relations for 
r{k). Typical dispersion relations are shown in Fig. 1. 
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FIG. 1: Instability growth rates as a function of wave number for 
^ = 10 and 6„ = 7l/8. Note that both growth rates vanish at ^: = and 
have a maximum T* --^ m/j / 10 at /:* --^ mo/3. The fact that they have 
a maximum means that at early times the system will be dominated 
by unstable modes with spatial frequency l/k* which grow at a 
rate T* . 



IV. DISCRETIZED HARD-LOOP DYNAMICS 

It is possible to go beyond an analysis of gluon polariza- 
tion tensor to a full effective field theory for the soft modes 
and then solve this numerically. The effective field theory for 
the soft modes that is generated by integrating out the hard 
plasma modes at one-loop order and in the approximation 
that the amplitudes of the soft gauge fields obey A ^ |p|/^ 
is that of the gauge-covariant collisionless Boltzmann-Vlasov 
equations [28]. In equilibrium, the corresponding (nonlo- 
cal) effective action is the so-called hard-thermal-loop effec- 
tive action which has a simple generalization to plasmas with 
anisotropic momentum distributions [12]. For the general 
non-equilibrium situation the resulting equations of motion 
are 

FMAy = [vD{A)]W^ix;y), (4) 

where / is a weighted sum of the quark and gluon distribution 
functions [12] and v^' = /5^/|p| = (1, v). 

These equations include all hard-loop resummed propaga- 
tors and vertices and are implicitly gauge covariant. At the 
expense of introducing a continuous set of auxiliary fields 
W|3(ji:;v) the effective field equations are also local. These 
equations of motion are then discretized in space-time and 
V, and solved numerically. The discretization in v-space cor- 
responds to including only a finite set of the auxiliary fields 
W|3(x;v,) with 1 < / < Nw- For details on the precise dis- 
cretizations used see Refs. [17, 18]. 
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FIG. 2: Plot of typical energy densities observed in non-expanding 
three-dimensional hard-loop simulation of the soft-fields generated 
in an anisotropic plasma with t, = 10. Shows transition from expo- 
nential growth with preference for transverse magnetic fields to linear 
isotropic growth. Inset shows same data on logarithmic vertical axis. 
To obtain physical units we have used a, = 0.3 and Qs = 1.5 GeV. 



A. Discussion of Hard-Loop Results 

During the process of instability growth the soft gauge 
fields get the energy for their growth from the hard particles 
and, of course, the total energy is conserved. In an Abelian 
plasma the energy deposited in soft fields grows exponen- 
tially until the energy in the soft fields is of the same order as 
the energy remaining in the hard particles at which point the 
back-reaction of the fields on the particle motion causes rapid 
isotropization. As mentioned above in a non- Abelian plasma 
the situation is quite different and one must rely on numeri- 
cal simulations due to the presence of strong gauge field self- 
interactions. 

In Fig. 2 I have plotted the time dependence of the chromo- 
magneticZ-electric energy densities obtained from a 3h-1 di- 
mensional from a typical HL simulation run initialized with 
"weak" random color noise with ^ = 10 [18]. The inset shows 
the data on a logarithmic scale. As can be seen from this fig- 
ure at f ~ 2.5 fm/c there is a change from exponential to lin- 
ear growth. Another interesting feature of the isotropic linear 
growth phase is that it exhibits a cascade of energy pumped 
into the unstable soft modes to higher energy plasmon like 
modes. This is demonstrated in Fig. 3 which shows the soft 
gauge field spectrum as a function of momentum at different 
simulation times in the saturated linear regime along with the 
estimated scaling coefficient of the spectrum [20]. 

Note, in addition, Ref. [20] showed that there were non- 
perturbatively large chaotic Chern-Simons number fluctua- 
tions in the linear growth phase. This should be contrasted to 
an Abelian theory or weak-field non- Abelian theory in which 
Chern-Simons number can fluctuate around zero but does not 
give large values when it is averaged over time. This is further 
indication of the non-perturbative physics associated with the 
chromo-Weibel instability. 
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FIG. 3: Field mode spectrum for SU(2) runs showing saturation of 
soft field growth at / ~ 1/g^ and an associated cascade of energy 
to the UV as the simulation time increases. Figure adapted from 
Ref. [20]. 



From Fig. 2 we can conclude that the chromo-Weibel in- 
stability will be less efficient at isotropizing a QCD plasma 
than the analogous Weibel instability seen in Abelian plasmas 
due to the slower than exponential growth at late times. On 
the positive side, from a theoretical perspective "saturation" 
at the soft scale implies that one can still apply the hard-loop 
effective theory self-consistently to understand the behavior 
of the system in the linear growth phase. 

One caveat is that the latest published HL simulation re- 
sults [17, 18] are for distributions with a finite o(l^lO) 
anisotropy due to computational limitations and saturation 
seems to imply that for weak anisotropics field instabilities 
will not rapidly isotropize the hard particles. This means, 
however, that due to the continued expansion of the sys- 
tem the anisotropy will increase. It is therefore important 
to understand the behavior of the system for more extreme 
anisotropics. More importantly it is necessary to study the 
hard-loop dynamics in an expanding system. Naively, one ex- 
pects this to change the growth from exp(T) to exp(-\/T) at 
short times but there is no clear expectation of what will hap- 
pen in the linear regime. A significant advance in this regard 
has occurred recently for a U(l) gauge theory [19]. Work is 
underway to do the same for non-Abelian gauge theories. 

V. WONG-YANG-MILLS EQUATIONS 

It is also to possible to go beyond the hard-loop approxima- 
tion and solve instead the full classical transport equations in 
three dimensions [23]. The Vlasov transport equation for hard 
gluons with non-Abelian color charge q" in the coUisionless 
approximation are [29, 30], 

f W - gq"F;^-d\ ~ 8UcA';,q%f]fix,p,q) = . (5) 

Here, f{t,x,'p,q^) denotes the single-particle phase space dis- 
tribution function. 



The Vlasov equation is coupled self-consistently to the 
Yang-Mills equation for the soft gluon fields, 

D,F^"'=r=g I ^dqqv'fit,x,p,q) , (6) 

with v" = ( 1 , p/p) . These equations reproduce the "hard ther- 
mal loop" effective action near equilibrium [24-26]. How- 
ever, the full classical transport theory (5,6) also reproduces 
some higher «-point vertices of the dimensionally reduced ef- 
fective action for static gluons [27] beyond the hard-loop ap- 
proximation. The back-reaction of the long-wavelength fields 
on the hard particles ("bending" of their trajectories) is, of 
course, taken into account, which is important for understand- 
ing particle dynamics in strong fields. 

Eq. (5) can be solved numerically by replacing the continu- 
ous single-particle distribution /(x,p,q) by a large number of 
test particles: 

/(x,p,q) = J_^8(x-Xi(t)) 

^Vtest / 

x(27:)35(p-pi(t))5(q''-q?(t)), (7) 

where Xi(t), Pi(t) and q"{t) are the position, momentum, and 
color phase-space coordinates of an individual test particle 
and Mest denotes the number of test-particles per physical par- 
ticle. The Ansatz (7) leads to Wong's equations [29, 30] 

^ - (E'' + VixB='), (9) 
dt 

^ = /§v^[A^,qi], (10) 

J"' = (11) 

^Vtest i 

for the /-th test particle. The time evolution of the Yang- 
Mills field can be followed by the standard Hamiltonian 
method [31] in A*^ = gauge. For details of the numerical 
implementation used see Ref. [23]. 

In Fig. 4 1 present the results of a three-dimensional Wong- 
Yang-Mills (WYM) simulation published in Ref. [23]. The 
figure shows the time evolution of the field energy densities 
for SU{2) gauge group resulting from a highly anisotropic 
initial particle momentum distribution. The behavior shown 
in Fig. 4 indicates that the results obtained from the hard- 
loop simulations and direct numerical solution of the WYM 
equations are qualitatively similar in that both show that for 
non-Abelian gauge theories there is a saturation of the energy 
transfeiTed to the soft modes by the gauge instability. Al- 
though I don't show it here the corresponding Coulomb gauge 
fixed field spectra show that the field saturation is accom- 
panied by an "avalanche" of energy transferred to soft field 
modes to higher frequency field modes with saturation occur- 
ring when the hardest lattice modes are filled [23]. 

A more thorough analytic understanding of this ultravio- 
let avalanche is lacking at this point in time although some 
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FIG. 4: Time evolution of the field energy densities for SU (2) gauge 
group resulting from a highly anisotropic initial particle momentum 
distribution. Simulation parameters are L = 5 fm, Phard = 16 GeV, 
g^ng = 10/fm^ m„ = 0.1 GeV. 



advances in this regard have been made recently [32]. Ad- 
ditionally, since within the numerical solution of the WYM 
equations the ultraviolet modes become populated rapidly this 
means that the effective theory which relies on a separation 



between hard (particle) and soft (field) scales breaks down. 
This should motivate research into numerical methods which 
can be used to "shuffle" field modes to particles when their 
momentum becomes too large and vice-versa for hard parti- 
cles. Hopefully, using such methods it will be possible to sim- 
ulate the non-equilibrium dynamics of anisotropic plasmas in 
a self-consistent numerical framework which treats particles 
and fields and the transmutation between these two types of 
degrees of freedom self-consistently. 



VI. OUTLOOK 

An important open question is whether quark-gluon plasma 
instabilities and/or the physics of anisotropic plasmas in gen- 
eral play an important phenomenological role at RHIC or 
LHC energies. In this regard the recent papers of Refs. [33- 
36] provide theoretical frameworks which can be used to cal- 
culate the impact of anisotropic momentum-space distribu- 
tions on observables such as jet shapes and the rapidity de- 
pendence of medium-produced photons. A concrete example 
of work in this direction is the recent calculation of photon 
production from an anisotropic QGP [37]. The results of that 
work suggest that it may be able possible to determine the 
time-dependent anisotropy of a QGP by measuring the rapid- 
ity dependence of high-energy medium photon production. 
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